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ABSTRACT 

We present an ultradeep A"s-band image that covers 0.5 x 0.5 deg 2 centered on the Great Observatories 
Origins Deep Survey-North (GOODS-N). The image reaches a 5 a depth of K s ,ab = 24.45 in the GOODS-N 
region, which is as deep as the GOODS-N Spitzer Infrared Array Camera (IRAC) 3.6 pm image. We present 
a new method of constructing IRAC catalogs that uses the higher spatial resolution K s image and catalog as 
priors and iteratively subtracts fluxes from the IRAC images to estimate the IRAC fluxes. Our iterative method 
is different from the \ 2 approach adopted by other groups. We verified our results using data taken in two 
different epochs of observations, as well as by comparing our colors with the colors of stars and with the colors 
derived from model spectral energy distributions (SEDs) of galaxies at various redshifts. We make available 
to the community our WIRCam A^-band image and catalog (94951 objects in 0.25 deg 2 ), the Interactive Data 
Language (IDL) pipeline used for reducing the WIRCam images, and our IRAC 3.6 to 8.0 pm catalog (16950 
objects in 0.06 deg 2 at 3.6 /im). With this improved Ks and IRAC catalog and a large spectroscopic sample 
from our previous work, we study the color-magnitude and color-color diagrams of galaxies. We compare 
the effectiveness of using Ks and IRAC colors to select active galactic nuclei (AGNs) and galaxies at various 
redshifts. We also study a color selection of z = 0.65-1.2 galaxies using the Ks, 3.6 /im, and 4.5 /im bands. 
Subject headings: catalogs — cosmology: observations — galaxies: evolution — galaxies: formation — galax- 
ies: high-redshift — infrared: galaxies 



1, INTRODUCTION 

Deep near-infrared (NIR) imaging is essential for under- 
standing galaxy evolution at high redshifts. One of its key 
applications is the determination of galactic mass. NIR 
luminosity is the best tracer of galactic stellar mass be- 
cause it is less affected by extinction and young stars. This 
role is traditio nally played by deep ground-based A"-band 
imaging (e.g.,[C owie et a lTi994tlKauffmann & Chariot! 1998t 
iBrinchmann & E llis 2000]). An important reason for this 
is that the A"-correction at A"-band is small across a broad 
range of redshifts and is relatively invariant against galaxy 
type s compared to the o ther two NIR bands (J and H, see, 
e.g., iKeenan et al.ll2010l) . Recently, the high sensitivity of 
the Spitzer Infrared Array Camera (IRAC, iFazioet al .1120041) 
adds additional powerful wavebands for this purpose (e.g., 
Fontan a et al.ll2006t iPerez-Gonzalez et al.ll2008t [Eisner et al. 



20081) . However, ground-based A"-band imaging remains im- 
portant because of the ease of observations. 

In order to study galaxy formation and evolution, we have 
been carrying out deep NIR imaging in the Great Observato- 
ries O rigins Deep Survey-North (GOODS-N iGiavalisco et ail 
l200l . In iBarger. Cowie. & Wangl (120081 hereafter B08) 
we presented a deep A^-band catalog in the GOODS-N ob- 
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tained with the Wide-field InfraRed Camera (WIRCam) on 
the 3.6 m Canada-France-H awaii Telescope (C FHT). The 
WIRCam image was used by Wang et al] d2007l) to study an 
optically faint submillim eter galaxy. The catalog was used by 
ICowie & Barged (12008b to construct a mass-selected sample 
of galaxies at z < 1.5 with an AB magnitude limit of 23.4. 
We have since added more data to the WIRCam K s imaging 
and slightly improved the reduction. In this paper we present 
the new image and a new Ks catalog. In order to achieve a 
more complete sampling of galaxy mass at high redshifts, we 
have also constructed a new 3.6 to 8.0 pm catalog based on 
Spitzer IRAC images in the GOODS-N (M. Dickinson et al., 
in preparation), which we also present here. 

Colors in the Spitzer IRAC bands are used to study 
the r edshifts of galaxies, especially optically faint ones 
(e.g.. iPope et alj|2006t IWang et al.ll2007t iDevlin eTaT]|2009b 
Mars den et al.ll2009l) . This is largely based on the rest-frame 
1 .6 /im bump in galactic spectral energy distributions (SEDs), 
which is cause d by the H~ opacity minimum in the stel - 
lar photosphere ( Simpson & Eisenhardt 1999; Sawicki 20021) . 
IRAC colors are also u sed for separating active galactic nu clei 
(A GNs) from galaxies dLacv et al.l2004tlStern et al.ll2005l se e 
also lHatziminaoglou et al.l2005USaiina. Lacv. & Scottfl 2005). 
In B08 we studied the effectiveness of using IRAC colors to 
select AGNs and galaxies at various redshifts based on a large 
spectroscopic sample, and we pointed out the limitations in 
these selections. With our improved color measurements in 
the IRAC and Ks bands, we revisit and discuss these issues. 

We describe the observations in § f2] the data reduction in 
§ [3] the reduction quality in § |4] and a K s selected catalog 
in § [5] We present in § [6] a new method of constructing 
IRAC catalogs based on the A^-band image and catalog, and 
we evaluate the performance of our new method. In § [7] we 
describe the general properties of the Ks and IRAC catalog, 
including its overlap with other multiwavelength catalogs, the 
Ks and IRAC color-magnitude diagrams and color-color dia- 
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grams. We summarize our results in §[8] The images and cat- 
alog data from this work will be available online. All fluxes 
in this paper are f v . All magnitudes are in AB system, un- 
less otherwise stated, where an AB magnitude is defined as 
AB = 23.9- 2.5 log(Jy). 

2. WIRCAM IMAGING OBSERVATIONS 

Ultradeep A^-band imaging observations of the Hubble 
Deep Field-North (HDF-N), including the GOODS-N and 
flanking fields, were carried out with WIRCam on the CFHT. 
WIRCam consists of four 2048 x 2048 HAWAII2-RG detec- 
tors covering a field of view of 20' x 20' with a 0"3 pixel 
scale. The observations were made in semesters 2006A, 
2007 A, and 2008A by two groups. A Canadian group led 
by Luc Simard obtained 11.5 hr of integration in 2006A. Our 
group based in Hawaii obtained 37.9 hr of integration between 
2006A and 2008A, bringing the total integration to 49.4 hr. To 
recover the gap between the sensors and bad pixels, all expo- 
sures were dithered with various offsets ranging from ±0.'5 to 
±1.'5 along R.A. and Dec. The dither patterns were changed 
after each semester in order to minimize artifacts caused by 
flat fielding or by sky subtraction. After each 0.5 to 1 hr of 
dithered observations, the dither centers were offset by up 
to ±6' to further broaden the area coverage and smooth out 
the distribution of integration time. The observations were 
carried out by the CFHT in queue mode with weather mon- 
itoring. Nearly all of the data were taken under photometric 
conditions. A small fraction (< 10%) of the data were taken 
under thin cirrus, so we have down-weighted them (see next 
section). Typical seeing for the observations is between 0."7 
and 0."8 (FWHM). Images taken under > 1" seeing were not 
used in the reduction. 

3. DATA REDUCTION 

In order to reduce imaging data taken with WIRCam and 
similar cameras, we have developed an Interactive Data Lan- 
guage (IDL) based reduction pipeline — SIMPLE Imaging and 
Mosaicking PipeLinE (SIMPLEfl Here we describe the SIM- 
PLE reduction of the WIRCam data. 

3.1. Removal of Instrumental Features 

The WIRCam images were reduced chip by chip. Because 
of the rapidly varying sky color in the NIR, only dithered 
images taken within ~ 0.5 hr were grouped and reduced to- 
gether. A group of raw images were first corrected for non- 
linearity in the instrumental response with correction coeffi- 
cients provided by the CFHT (C.-H. Yan, 2008, priv. comm.). 
The images were then normalized, median-combined to form 
a sky flat field, and flattened with that. In the flattened im- 
ages, objects were detected and masked back in the raw im- 
ages. The masked raw images were normalized again by tak- 
ing into account the locations of the masked pixels and then 
median-combined to form a second sky flat. This was used as 
the final flat field for flattening the raw images again. Because 
there was still some variation in sky color within the ~ 0.5 hr 
period, the sky flat did not always perfectly flatten the images. 
To remove the residual gradients in the flattened images, ob- 
jects were detected again and masked. A 5th-degree polyno- 
mial surface was fitted to each masked image and subtracted 
from the unmasked image. 

7 The SIMPLE package is available at 
http://www.asiaa.sinica.edu.tw/~whwang/idl/SIMPLE 



On each WIRCam detector, there is crosstalk between each 
of the 32 readout channels (2048 x 64 pixels each), as well 
as within each of the four video boards (8 channels each). 
For every image, SIMPLE median-combined the 32 2048 x 
64 stripes and subtracted the combined stripe from each of 
the original ones to remove the 32-channel crosstalk features. 
The 8-channel crosstalk features are relatively weak and only 
appear around bright objects after combining several hours of 
images. SIMPLE repeated the median combination around 
only bright objects and only within the affected 8 channels, 
and then subtracted the median image to remove the 8-channel 
features. All the median combinations here were performed 
after detected objects were masked. 

The above procedures cleanly removed most of the 32 and 
8-channel crosstalk. Very low level 8-channel crosstalk fea- 
tures are still visible around the brightest ~ 5 stars in our 
field-of-view after combining the 39 hr of data taken in 2006 
and 2007. This is probably caused by the coupling between 
crosstalk, flat field pattern, and sky background structures in 
the images. There is no trivial way to solve this, and we 
did not attempt to further remove this very low level effect. 
Thanks to the WIRCam instrument team, the crosstalk in 2008 
had improved behavior and no crosstalk residual effect is ob- 
served after the above 32 and 8-channel removal procedure. 

3.2. Correction for Distortion and Astrometry 

The flattened, sky subtracted, and crosstalk removed im- 
ages were then used for deriving the pointing offsets between 
the dithered images and the distortion function of the optics. 
In order to obtain accurate offsets, SIMPLE uses the package 
SExtractor (Bertin & Arnoutslll996l) to detect objects and to 
determine their pixel coordinates. Then offsets are derived by 
averaging the displacements of the same objects that appear 
in different images. Because of the optical distortion, objects 
in different parts of the images have slightly different dither 
displacements. This differential displacement is used by SIM- 
PLE to deriv e the distortion function, f ollowing the technique 
developed by I Anderson & King ( 2003|). 

SIMPLE uses cubic polynomials to fit the distortion func- 
tions, i.e., 



x' =F(x,y) = ^2f u x'y J , 
y' = G(x,y) = ^ gi jx'y ] , 



(1) 

(2) 



where x and y are the undistorted coordinates of objects, and x' 
and y' are the distorted coordinates, F and G are the distortion 
functions, and i + j <3. To first order, the displacements of 
objects between dithered images can then be expressed as the 
expansions of the distortion functions: 



dx 



dy 



dx 



dy 



(3) 
(4) 



where Ax' and Ay' are the measured dither displacements of 
objects, which are functions of positions, and Ax and Ay are 
the pointing offsets of the dithering. 

Under the above approximation, the measured displace- 
ments of objects are the first order derivatives of the distortion 
functions. For n objects in m dithered images, there are ap- 
proximately n x (m- 1) sets of linear equations for each of Ax' 
and Ay' . The systems of linear equations were solved for the 
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coefficients of F and G with a least-squares method. Initially, 
the mean values of Ax' and Ay' were used as rough estimates 
of Ax and Ay. After F and G were integrated, they were used 
to correct for x 1 and y' to derive better estimates of the dither 
offsets Ax and Ay. Then improved coefficients of F and G 
were solved with the new Ax and Ay. For WIRCam, one 
such iteration seems to provide sufficiently good distortion 
functions. Since the optical distortion was derived within a 
~ 0.5 hr sequence of dithered images, the above approach has 
the advantage of accounting for the long-term time-dependent 
flexure of the telescope. It also does not require any external 
information and thus can easily be applied to any target field. 

To determine the absolute astrometry, we used several dif- 
ferent catalogs in the HDF-N region. We first combined the 
catalogs of the Sloan Digital Sky SurvejQ (SDSS) Data Re- 
lease 6 around the HDF-N region and the vl. O GOODS-N Ad- 
vance d Camera for Surveys (ACS) catalog (Gia valisco et al.l 
2004), after correcting for the well-known 0."38 offset be- 
tween the GOODS-N ACS frame and the radio frame. The 
combined catalog was cross-checked with the source posi- 
tions in the Ve r y Lar ge Array (VLA) 1.4 GHz catalog of 
Biggs & Ivisonl (120061) to ensure consistency with the radio 
frame. The comparison between the ACS and the VLA as- 
trometry is only good to several tens of milli-arcseconds, be- 
cause of the differences in optical and radio morphologies of 
distant galaxies. This sets the limit for our absolute astrome- 
try. 

We then reduced several hours of the best quality WIRCam 
A^-band data and matched the WIRCam object positions to 
the above SDSS/ACS astrometric catalog. We added bright 
and compact WIRCam K s objects to the SDSS/ACS catalog 
to form a dense astrometric catalog for reducing the full set 
of WIRCam Ks images. The final catalog has ~ 8700 ob- 
jects over a diameter of 1.3 degree, approximately 5700 of 
which are concentrated in the central 0.6 deg 2 region. In all 
the reductions, when individual images were warped to cor- 
rect for distortion and projected onto a common sky plane, 
the objects detected in the images were forced to match the 
positions in the astrometric catalog. The registration of the 
images onto the catalog was performed at the sub-pixel level 
with a bilinear resampling. There was only one resampling for 
each exposure in the entire reduction, which was a fifth degree 
polynomial including both the distortion correction and a tan- 
gential sky projection. This one-step resampling minimized 
the image smearing. 

3.3. Mosaicking and Photometric Calibration 

8 Funding for the SDSS and SDSS-II has been provided by the Alfred P. 
Sloan Foundation, the Participating Institutions, the National Science Foun- 
dation, the U.S. Department of Energy, the National Aeronautics and Space 
Administration, the Japanese Monbukagakusho, the Max Planck Society, and 
the Higher Education Funding Council for England. The SDSS Web Site is 
http://www.sdss.org/ 

The SDSS is managed by the Astrophysical Research Consortium for the 
Participating Institutions. The Participating Institutions are the American 
Museum of Natural History, Astrophysical Institute Potsdam, University of 
Basel, University of Cambridge, Case Western Reserve University, Univer- 
sity of Chicago, Drexel University, Fermilab, the Institute for Advanced 
Study, the Japan Participation Group, Johns Hopkins University, the Joint 
Institute for Nuclear Astrophysics, the Kavli Institute for Particle Astro- 
physics and Cosmology, the Korean Scientist Group, the Chinese Academy 
of Sciences (LAMOST), Los Alamos National Laboratory, the Max-Planck- 
Institute for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics 
(MPA), New Mexico State University, Ohio State University, University of 
Pittsburgh, University of Portsmouth, Princeton University, the United States 
Naval Observatory, and the University of Washington. 




filter and the dashed curve is the 2MASS Ks filter. 

Before the distortion-corrected images were combined to 
form a deep mosaic, the images were calibrated for extinction 
variation, cosmic rays were removed, and the images were 
weighted. The relative extinction variation was corrected us- 
ing fluxes measured by SExtractor on high S/N objects in the 
images. Then pixels in different images that have the same 
sky coordinates were compared for rejecting cosmic rays with 
sigma clipping. The clipping rate was kept to be approxi- 
mately 1 % so only an insignificant amount of useful data were 
discarded. The sigma clipping allows weighted mean combi- 
nations of the images, which produces 25% higher S/N com- 
pared to median combinations when all images have equal 
weights. The images were optimally weighted according to 
the atmospheric transparency (?7 a t m , derived from the flux cal- 
ibration, see below), integration time (r exp ), and background 
brightness (BG). Furthermore, each pixel was weighted by 
its relative efficiency (Qe, proportional to flat field). Thus the 
final weight applied to a pixel is: 

w(x,y) oc ?7 2 tm 7; X pBG" 1 2E(^,3 ; )- (5) 

Since objects in the images have a wide range of sizes, we 
did not weight the images with seeing. The relative extinction 
corrected, cosmic ray removed, and weighted images were 
then mean-combined to form a deep mosaic. 

On the mosaicked images, object fluxes were measured 
again with SExtractor using 5" diameter apertures, and the 
results were cross- correlated with the T wo Micron All Sky 
surve >0 (2MASS, Skrutskie et al. 2006) point-source cata- 
log. We used 2MASS point sources with K s = 14.44-16.34 
to calibrate the absolute flux of the WIRCam image. Objects 
brighter than this range are affected by WIRCam nonlinearity 
(despite the correction for linearity) and objects fainter than 
this range show a significant selection effect in the 2MASS 
fluxes. We adopted the 2MASS "default magnitudes," which 
attempt to accou nt for the total fluxes of the point sources 
dCutri et alj|2003l) . Because the 2MASS K s and WIRCam K s 
filter transmission curves are very similar (Figure Q]), we di- 
rectly adopted the 2MASS Ks magnitudes and did not attempt 
to account for the slight difference between the two systems. 
We adopted the 666.7 Jy Ks Vega zero-magnitu de provided 
by 2MASS dCohen, Wheaton, & Megeafhll2003l) . Calibrated 
mosaic images from all the dither sets and from the four chips 
were combined into a final wide-field, ultradeep image. All 
the weights described above were carried into the final image 
in the form of effective integration time. We found that de- 
spite the large variation in sky background brightness, noise 
in the final image scales nicely with the inverse square-root 
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FIG. 2. — Survey field layout. Grayscale shows the effective integration 
time of our WIRCam K s imaging. White contours show 80%, 60%, 40% 
and 20% of the maximum effective integration time. Black solid polygons 
show the GOODS-N ACS region and the HDF-N-proper. Black dashed con- 
tour shows the GOODS-N IRAC region. Black rectangle shows the Chandra 
Deep Field-North region. 
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FIG. 3. — Cumulative area in the WIRCam Ks image vs. effective integra- 
tion time. The maximum effective integration time in the map is normalized 
to be 1 .0. The effective integration time accounts for variation in efficiency of 
pixels, atmospheric extinction, and background brightness (see text). In shal- 
low and intermediate depth regions where confusion is less important, noise 
very closely scales as inverse square-root of the effective integration time. 



of the effective integration time in intermediate and shallow 
regions where confusion noise is less important. This means 
that our weighting scheme provides excellent S/N optimiza- 
tion. 

Figure |2] shows the integration time distribution and layout 
of our survey field. Figure [3] shows the cumulative area in 
the final mosaic as a function of effective integration time. 
The final combined mosaic covers a square area of 35' x 35', 
centered at a = 189°. 227, 6 = 62°. 239 (J2000.0). The central 
~ 18' region has more than 20 hr of integration at each posi- 
tion and fully covers the GOODS-N ACS area. The image has 
uniform image quality across the entire field, with a FWHM 
of~0"8. 

9 2MASS is a joint project of the University of Massachusetts and the 
Infrared Processing and Analysis Center/California Institute of Technology, 
funded by the National Aeronautics and Space Administration and the Na- 
tional Science Foundation. 
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FIG. 4. — Relative astrometric offsets between the WIRCam Ks image and 
the SDSS and GOODS-N ACS catalogs. Approximately 4000 well-detected 
compact sources (S/N > 20, FWHM < If 2) are plotted. Squares show the 
offsets of the WIRCam positions with respect to the SDSS positions, and dots 
show those with respect to the ACS positions. 

4. REDUCTION QUALITY 
4.1. Astrometry 

The SIMPLE reduction of the WIRCam Ks image forced 
the positions of detected objects to match those in the astro- 
metric catalog formed with the SDSS and the GOODS-N ACS 
catalogs. The resultant astrometry is excellently constrained 
by the ACS catalog in the central region and reasonably con- 
strained by the SDSS catalog in the outer region. To quantify 
the (relative) astrometric accuracy, we selected well-detected 
compact sources (S/N > 20 and FWHM < 1 ."2) from the SEx- 
tractor catalog (Section 01. Nearly 4000 such sources have 
counterparts in the SDSS or GOODS-N ACS catalog. Fig- 
ure|4]and Figure|5]show the offsets between the WIRCam and 
the SDSS/ACS positions. There are no observable systematic 
offsets between the WIRCam positions and both the SDSS 
and ACS catalogs. However, the astrometry is clearly better 
in the central part of the image where ACS sources were used 
to constrain the astrometric solution. The rms dispersions of 
the offsets are 0."05 between WIRCam and ACS along both 
the R.A. and Dec. axes on nearly 4000 sources. The rms off- 
sets between WIRCam and SDSS are larger, 0."08 for RA and 
0"10 for Dec, on nearly 460 sources. We found that the SDSS 
and ACS catalogs have similar ~ 0."1 rms offsets between 
each other. This suggests the larger observed offsets between 
the WIRCam and SDSS positions may primarily be a conse- 
quence of an uncertainty internal to the SDSS catalog. 

The above measured rms offsets include both systematic 
and random errors. For example, the measured 0."05 rms 
offset between WIRCam and ACS is slightly larger than 
the nominal positional uncertainty for FWHM ~ 0"8 and 
S/N > 20, which should be < 0."04. To estimate the system- 
atic error in astrometry, we repeated the same measurements 
with higher S/N cuts and lower FWHM cuts. The rms offset 
decreases initially but stops decreasing when it approaches 
0"03, even when the selected sources have extremely high 
S/N (e.g., > 200) and have FWHM comparable to the seeing. 
We conclude that there is an ~ 0."03 systematic uncertainty in 
our WIRCam astrometry in the GOODS-N ACS region. Such 
an error includes both the internal error in the ACS catalog 
and the error introduced by the image registration of SIM- 
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FIG. 5. — Same as Figure|4] but as functions of RA and Dec. 

PLE. On the other hand, the rms offset between WIRCam and 
SDSS is — 0."1. Since all SDSS sources are detected in the 
WIRCam image with S/N > 100 (median S/N is 2300), ran- 
dom noise is negligible here. We conclude that the systematic 
error is ~ 0."1 in outer regions where ACS positions are un- 
available. 

4.2. Photometry 

The absolute flux scale of our WIRCam K s image was cal- 
ibrated with 2MASS point sources in the image. Figure [6] 
shows the comparison between the 2MASS fluxes and the cal- 
ibrated WIRCam fluxes derived with 5" diameter apertures. 
The ratio between the two is only flat in a narrow flux range. 
We chose K^ab = 14.44-16.34 (the two vertical lines in Fig- 
ure [6]) for the calibration of the WIRCam image. Outside 
this range, there are nonlinearity problems with the WIRCam 
sensors (bright end) and selection effects in 2MASS (faint 
end). In the chosen range for calibration, the error- weighted 
mean of the flux ratios (^2MASs/-FwiRCam) for 70 sources is 
1.0005 ±0.5%, indicating that the overall absolute calibra- 
tion of the WIRCam Kg image is good to ~ 0.5%. We note 
that the error bars in Figure|6]are dominated by 2MASS flux 
errors. Thus, the calibration of our WIRCam image is lim- 
ited by the photometric quality of 2MASS. It is not possible 
to achieve better results unless we make frequent monitoring 
of standard stars during our observations, which is quite ex- 
pensive in terms of observing time. 

We also investigated the photometric uniformity in our 
WIRCam Ks image. The above value of 0.5% is the mean cal- 
ibration error over the entire image. It is possible that the cali- 
bration of certain smaller regions is worse than this value. We 
subdivided the WIRCam image into four quadrants, roughly 
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FIG. 6. — Comparison of the 2MASS Kg fluxes with the calibrated 
WIRCam Kg fluxes. The noise in the WIRCam fluxes is plotted as horizontal 
error bars but is too small to be visible for the entire flux range and therefore 
makes a negligible contribution to the vertical error bars. This implies a cali- 
bration limited by the photometric quality of 2MASS. Vertical lines show the 
flux range used for the calibration. Sources outside the two lines are not used 
because of nonlinearity in WIRCam fluxes and selection effects in 2MASS. 
Within the two vertical lines, data that are > 3cr away from the mean are also 
excluded in the calibration. 

corresponding to data obtained from the four detector chips, 
and then compared the WIRCam fluxes with the 2MASS 
fluxes as we did for the entire image. For quadrants 1 
through 4, the error-weighted flux ratios are 1.010 ±0.9% 
(16 sources), -0.985 ±0.9% (19 sources), 1.009 ± 1.3% (16 
sources), and 1.011 ±0.9% (22 sources), respectively. In all 
of these cases, based on the numbers of available 2MASS 
sources and their 2MASS flux errors, the expected calibra- 
tion errors are roughly 1%. The standard deviation of the four 
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TABLE 1 
SExtractor Parameters 



Parameter Value 



DETECT_MINAREA 




2 


DETECT_THRESH 




1.25 


ANALYSIS_THRESH 




1.25 


FILTER 




Y 


FILTER NAME 


gauss_1.5 


3x3.conv 


DEBLEND NTHRESH 




64 


DEBLEND MINCONT 




0.00001 


CLEAN 




Y 


CLEAN PARAM 




0.1 


SEEING FWHM 




0.8 


BACK SIZE 




24 


BACK FILTERSIZE 




6 


BACK TYPE 




AUTO 


BACKPHOTO TYPE 




LOCAL 


BACKPHOTO THICK 




24 


WEIGHT TYPE 


MAP_ 


WEIGHT 


WEIGHT THRESH 
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measured offsets is 1.3%, which compares well with the ex- 
pected calibration error from 2MASS. This suggests that there 
is not an observable photometric gradient over size scales of 
~ 15'. The absolute fluxes of objects are still limited by the 
photometric quality of 2MASS at this scale. The measured 
1.3% variation is indeed excellent when compared to other K- 
band extragalactic deep sur veys over similar si ze scales (e.g., 
2% in the COSMOS survey. ICapak et alj|2007l) . 

5. K$ SOURCE CATALOG 

We used SExtractor (v2.5.0) to produce a Ks source catalog. 
The most important source extraction parameters in SExtrac- 
tor are summarized in TableQ] We provide our raw SExtractor 
Ks catalog in this data release. We caution that the detection 
parameters (the first three in the table) that we adopted are 
somewhat aggressive in detecting fai nt so urces, which may 
lead to more spurious sources (see § |5.2| i. The users of our 
data product may want to regenerate catalogs with more con- 
servative values, based on the nature of their work. 

Based on the raw catalog, we compiled a K$ selected cata- 
log, which is the primary product of this paper. We first es- 
timated the true noise level in the WIRCam K s image. We 
masked the detected objects and convolved the image with 
circular kernels of amplitude 1.0 (i.e., flux apertures). The 
local (within ~ 14") rms fluctuations in the convolved image 
were considered as errors in aperture photometry that include 
all possible sources of noise: photon noise, correlated noise 
between pixels, confusion noise from faint, undetected ob- 
jects, and uncertainty in the estimate of the small-scale back- 
ground. We compared the noise measured in this way with the 
photometric errors derived by SExtractor. We found that for 
all aperture sizes between l."5 and 5" in diameter, the ratio 
between the two errors is fairly insensitive to the aperture size 
and has a value of 1.28. Thus, we multiplied all of the flux 
errors in the SExtractor catalog by this factor to approximate 
the total noise. This does not affect the flux measurements. 

The next step is to estimate aperture corrections. In our data 
reduction we calibrated the flux on 2MASS point sources us- 
ing fixed 5" diameter apertures. Thus, SExtractor fluxes mea- 
sured with 5" apertures do not need any aperture corrections 
and we consider 5" aperture fluxes to be the best estimate 
of the total fluxes of the sources. This clearly does not ap- 
ply to the several most extended objects in the image, since 
their aperture corrections should be much larger than that for 
smaller objects (i.e., 2MASS point sources used in the cali- 



bration). Since they are not the major targets of our work, we 
did not attempt to obtain highly accurate fluxes for these very 
few sources. 

Despite that 5" aperture fluxes are the best estimate of total 
fluxes, their S/N is relatively low. We thus looked for suit- 
able fluxes measured with smaller apertures, and corrected the 
small aperture fluxes to 5" when it is necessary. We compared 
all other SExtractor flux options with the fixed 5" diameter 
aperture fluxes and found that the SExtractor "auto-aperture" 
fluxes agree excellently with them at all flux levels. In ad- 
dition, the auto-aperture fluxes have the best optimized S/N 
for sources with various morphologies and fluxes. We thus 
decided to adopt directly the SExtractor auto-aperture fluxes 
and errors in our final catalog. We did not apply any additional 
corrections except for the aforementioned factor of 1 .28 to the 
errors. 

The adoption of SExtractor auto-aperture fluxes works ex- 
cellently for nearly all sources. However, a very small num- 
ber of sources (0.7% of total) are undetected in auto-apertures 
but detected in fixed apertures. For such sources, we adopted 
their 3" fluxes and corrected the fluxes to 5" with a median 
correction factor from all other detected sources. In addition, 
after all the above corrections to the fluxes and flux errors, 
some sources in the SExtractor raw catalog are still unde- 
tected. Given that our error estimate is realistic and the flux is 
well optimized, we considered sources with S/N < 1 as unde- 
tected ( 1 .9% of total) and only included their flux errors in the 
final catalog for upper limits. To keep the source number and 
order in our final Ks catalog identical to those in the SExtrac- 
tor raw catalog, we did not remove such undetected sources. 
This is for easy comparisons with the SExtractor catalog, in 
case the users of our data need to carry out more consistency 
or quality check. 

We next evaluated the depth and completeness of our Ks 
WIRCam catalog. In the central 30' x 30' of the image where 
the sensitivity distribution is more uniform, the median lcr 
error among 3a sources is 0.15 /xJy, corresponding to a 3a 
limiting magnitude of 24.77. In the GOODS-N region, where 
the image is the deepest, the median la error is 0.12 pjy, 
corresponding to a 3a limiting magnitude of 25. 10. This is as 
deep as the IRAC 3.6 pm image in the same region. 

5.1. Depth and Completeness 

We estimated the completeness of our Ks catalog with 
Monte Carlo simulations. We note that completeness is a 
complex function of flux, morphology, crowdedness, and 
source extraction. The results here are meant only to provide 
a guide. Exactly how one should carry out the completeness 
computation depends on the purpose of the study, and we shall 
leave this to the users of our data. 

Here, we first randomly selected well-detected (S/N > 20) 
sources from the SExtractor catalog. The selected sources 
span a wide range of flux, size, and shape. We extracted these 
sources from the WIRCam image and artificially dimmed or 
brightened them to the flux range of interest. Then we ran- 
domly placed them in the WIRCam image and ran SExtractor 
to detect them. Each time we only placed a small number 
of sources in the image so the overall source density did not 
change by > 1%. A total of ~ 50,000 sources were randomly 
placed in the image. The SExtractor catalogs for simulated 
images were processed identically to the real catalog. We con- 
sidered a source recovered if it is detected at > 2a. We divided 
sources into four groups according to their half-light radii (/?/,) 
and calculated the recovery rate in each 0. 1 dex flux bin and 
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FIG. 7. — Completeness of the Kg WIRCam catalog vs. flux and half-light 
radius in the central 30' X 30' region (A) and in the GOODS-N IRAC region 
(B). Percentages given in the parentheses are fractions of objects in the four 
half-light radius groups in the K$ catalog. 

Ri, group. The results in the central 30' x 30' region are shown 
in Figure [7^. The 90% completeness limits for the two com- 
pact groups (R h ~ 0."47 and 0."54, comprising > 50% of the 
total sample) are 1.10 /iJy (7 3a) and 1.32 /iJy (8.8cr), respec- 
tively, equivalent to 23.80 and 23.60 magnitude. For the entire 
sample the mean 90% completeness weighted by the size dis- 
tribution is 1.38 /iJy (9.2cr, which is equivalent to 23.55 mag- 
nitude. If we restrict this to the deeper 0.06 deg 2 GOODS-N 
IRAC region (Figure |7p), then the mean 90% completeness 
limit becomes 0.94 /iJy (7.8cr), or 23.96 magnitude. 

The above results show that the present WIRCam imag- 
ing is the deepest in the A'-band over a substantial area. For 
example, the Very Large Telescope (VLT) NIR imaging in 
the GOODS-S is 90% complete at K s = 23.8 in a small area 
of < 10 arcmin 2 dGrazian et al.l 120061). On the other hand, 
our team (IWang. Bareer. & Cowid 2009) and Japanese groups 
(Kaiisawa et al. 2006J) have carried out extremely deep K s 
band imaging in the GOODS-N at the 8.2 m Subaru Tele- 
scope with a smaller 2048 x 2048 x 2 imager. Kaiisa wa et al.l 
(2009) combined all archived data to make a deep Ks image 
of the GOODS-N. They achieved 3a limiting magnitudes of 
25.5 and 26.5 in 103 arcmin 2 and 28 arcmin 2 , respectively. 
This Subaru image is deeper than our WIRCam image (3a of 
25.1 magnitude in the ~ 200 arcmin 2 GOODS-N region), but 
the area coverage is much smaller. 

5.2. Spurious Sources 

We took two approaches to estimate the spurious fraction, 
one based on an inverted image and the other based on multi- 
wavelength counterpart identification. First, we constructed 
an image with inverted background, and then ran through 
the same source detection and catalog compilation described 



TABLE 2 

Properties of the WIRCam Images 



Properties 



Ks Image 



Semesters 
P.I. 

Total Integration 
Field Center (J2000.0) 

R.A. 

Dec. 

Area (T exp > 2 hi) 
Seeing FWHM 
Calibration Uncertainty 
Astrometric Uncertainty 
90% Completeness Limit 
Spurious Fraction 



2006A, 2007A, 2008A 
L. Cowie, L. Simard 
49.4 hi- 

189° .227 
62° .239 
0.25 deg 2 
0."8 
1.3% 
0"03-0"1 
1.38 /jjy (23.55 mag.) 
< 3.2% 



above. This infers a spurious fraction of ~ 20% on sources 
with S/N > 3, which is unreasonably high. After examin- 
ing the image and the "sources," we believe this is because 
the crosstalk of WIRCam produces negative holes near bright 
objects, making the noise non-Gaussian and greatly overesti- 
mating the spurious fraction. One evidence for this is that, in 
relatively crosstalk-free regions, the spurious fraction inferred 
from the inverted image dramatically decreases to 5%. This 
value of 5% is an upper limit since there are still low-level 
crosstalk effects. 

We then looked for K s sources that are not detected in other 
images. In the region where the GOODS-N ACS and IRAC 
images overlap, there are 15112 Ks sources detected at 3cr. 
Among them, 489 (3.2%) sources do not have counterparts 
in the ACS catalog and are not detected at 3.6 and 4.5 ^tm at 
3 a in the IRAC catalog described in the next section. (See 
§ 17. II for a slightly different comparison.) The fraction 3.2% 
is also an upper limit since we cannot rule out real Ks sources 
being not detected in the ACS and IRAC images, especially 
given that the IRAC images are not substantially deeper than 
the Ks image. This 3.2% upper limit is tighter than the 5% 
limit set by the inverted image, consistent with the above sug- 
gested non-Gaussian noise. Because of this limitation, it is 
difficult to further constrain the spurious fraction of our K$ 
catalog. We suggest the users of our data to carry out inde- 
pendent source extraction if more conservative extractions are 
desired. 

In Table [2] we summarize the basic properties of our 
WIRCam Ks image described in this and previous sections. 

6. COLORS IN THE WIRCAM AND IRAC BANDS 

We next used the WIRCam K$ image and catalog as 
priors to analyze the fluxes of galaxies in the GOODS-N 
Spitzer IRAC images at 3.6-8.0 /im (M. Dickinson et al., 
in preparation). Sim ilar work has been done previously. In 
both the GOODS-N dLaidler et al.ll2007h and the GOODS-S 
dGrazian et alJ20 06). ACS z'-band images were used as priors 
to determine fluxes in the IRAC bands. While our Ks image 
has a lower resolution than the ACS z' images, the waveband 
is much closer to the IRAC bands, making it a better choice 
for this type of analysis. 

The primary goal is to better measure the colors of objects 
between the IRAC and Ks bands. The nearly 2" resolution of 
IRAC is significantly worse than the typically sub-arcsecond 
resolution in the optical and NIR, but it is not much larger 
than the sizes of galaxies. This makes color measurements 
extremely uncertain. Normally one would need to sacrifice 
resolution in the optical and NIR bands to achieve this. The 
approach we describe below does not have this disadvantage. 
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A secondary goal is to better separate IRAC fluxes of objects 
that are close to each other. There are nearly 20,000 K$ de- 
tected source in the GOODS-N IRAC area, but SExtractor 
can only directly detect 11,000 sources in the 3.6 /mi IRAC 
image. One reason for the smaller number of sources in the 
IRAC image is that many objects are blended under the IRAC 
resolution. We will demonstrate that our method provides rea- 
sonable estimates of the IRAC fluxes of lightly blended ob- 
jects. 

6.1. Methodology 

The basic idea is to convolve the K$ image of each galaxy 
individually with a scaled convolution kernel to construct a 
model image for an IRAC image. The scaling of the ker- 
nel varies between galaxies because galaxies have different 
IRAC-to-Aj flux ratios (i.e., colors). The set of scaling factors 
that produce the best model image are then the best estimates 
of the flux ratios. 

It is immediately obvious that the results would be sensi- 
tive to the evaluation of whether a model best approximates 
the IRAC image, as well as wheth er the chosen c onvo lution 
kernel well represen ts the PSF. In lLaidler et alj (12007b and 
Grazian et al. (2006), the model that produced the minimum 
X 2 was considered as the best model. In such analyses, x 2 is 
very sensitive to the mismatch of the PSF, and unfortunately, 
there is a huge uncertainty in the PSFs in the IRAC bands. 
This was pointed out by Grazian et al. and is also clearly 
shown in Figure [8] It is thus unclear to us whether or not x 2 
analyses are the best approach. Here we present an alternative 
approach based on the widely used imag e deconvo lution tech- 
nique in radio astronomy — CLEAN ( Ho gbo ml 19741) . Instead 
of minimizing the x 2 of the pixels, CLEAN determines the 
scaling of the kernels based on the fluxes of the image com- 
ponents. This is very attractive to us since our true goal is 
to be able to account for the fluxes detected in the IRAC im- 
ages. As long as the fluxes in an IRAC image can be cleanly 
removed by subtracting a model, we consider the model to be 
the best. Below we describe the basic steps in our analysis of 
the galaxy colors. We hereafter refer to our technique as RE- 
ALCLEAN, since it is done in real space instead of Fourier 
space. 

6.1.1. PSF Construction 

The GOODS-N IRAC images have pixel scales of 0."6, 
so we background subtracted and resampled them to the 
WIRCam 0."3 pixel scale. During the resampling, the IRAC 
images were also slightly warped such that the positions of 
objects precisely match those in the WIRCam imag43 To 
construct PSFs for the IRAC and WIRCam images, we se- 
lected ~ 20 bright but unsaturated, point-like, and isolated 
objects in these images. The images of these objects were 
aligned, stacked, and normalized to form the PSFs. The IRAC 
PSFs were then deconvolved by the WIRCam PSF with a 
maximum likelihood deconvolution of a few iterations. The 
deconvolved IRAC PSFs were used as the convolution kernels 
to convolve the WIRCam image. One may worry about the 
uncertainty in the deconvolution of the PSFs. We also tried an 
alternative approach that does not require the deconvolution 
of the IRAC PSF, in which we measured colors in the IRAC 
images that are convolved with the WIRCam PSF and i n th e 
Ks image that are convolved with the IRAC PSFs (see § 16.4b . 
We did not see a systematic difference between the two ap- 
proaches. We thus adopted the deconvolved IRAC PSFs as 



the convolution kernels. 

6.1.2. Object Identification 

In the background-subtracted WIRCam Ks image, for each 
detected object we identified pixels associated with that ob- 
ject. This was done with the "segmentation map" provided by 
SExtractor. This allowed us to apply different scaling factors 
to the convolution kernels on different objects. A minor draw- 
back here is that very faint pixels in the outer parts of galaxies 
would be missed. However, as long as these faint pixels do 
not have dramatically different colors from the brighter cores 
of the galaxies, our measured colors can still be considered 
as good estimates of the flux-weighted mean colors of galax- 
ies. Another minor drawback is pixels that contain emission 
from more than one object. This would only become prob- 
lematic when objects are not well resolved in the WIRCam 
image. Since this is the fundamental limit of the WIRCam 
observations, we did not attempt to resolve this issue. 

6.1.3. REALCLEANing Bright Objects 

To obtain model images that best approximate the observed 
IRAC images, we use our REALCLEAN method, which is es- 
pecially useful when determining fluxes of blended galaxies. 
In a REALCLEAN cycle, we first identify the brightest galaxy 
in the IRAC image. Then we remove a small fraction (25% 
at 3.6 /im down to 10% at 8.0 /im; this factor is often called 
the "gain") of its flux from the IRAC image by subtracting a 
model image constructed by convolving the WIRCam image 
with the above mentioned kernel. The flux is determined by 
placing a circular aperture at the object. We used aperture di- 
ameters of 4", 4", 5", and 6" for the IRAC 3.6, 4.5, 5.8 and 
8.0 /im bands, respectively. Despite the fact that such aperture 
photometry is affected by nearby objects, because we only 
look at the brightest object and because the gain is <C 1, we 
guarantee that we do not over-subtract the object. To speed up 
the REALCLEAN process, we subtract fractional fluxes from 
many objects in one REALCLEAN cycle. However, in each 
REALCLEAN cycle we only subtract objects that are suffi- 
ciently far away from each other that locally each subtracted 
object is the brightest. In the subsequent REALCLEAN cycle 
we look for the next brightest objects in the subtracted image 
and subtract small fractions of fluxes from them. We repeat 
this iterative procedure until no objects in the residual IRAC 
image have fluxes brighter than a few a. 

It is interesting to note that the above procedure is similar to 
the CLEAN deconvolution of radio interferom etric data (see a 
summary in Cornwell, Braun, & Briggs 1999) in many ways. 
To limit the degrees of freedom, CLEAN is often limited to 
CLEAN windows or boxes within which emission is known 
or guessed to be present. CLEAN also sometimes uses priors 
to improve the speed of convergence of the solution. These 
are exactly what our procedures do: we only subtract IRAC 
fluxes at the locations of Ks detected galaxies, and we use the 
WIRCam Ks images of the galaxies as the model. 

6. 1 .4. Final Residual Flux Minimization 

Once REALCLEAN reaches a few a, all objects have sim- 
ilar brightness and REALCLEAN is no longer efficient and 
necessary. Thus, we measured the residual IRAC fluxes of all 
of the WIRCam objects and subtracted them from the IRAC 
image all at once. This global subtraction unavoidably over- 
subtracts blended objects, especially those whose separations 

10 This is achieved with subroutines provided in the SIMPLE package. 
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are less than the flux aperture sizes. Since the global subtrac- 
tion is carried out only after the objects are REALCLEANed 
to a few a, the amount of over-subtraction is small. Neverthe- 
less, the over-subtraction can be further minimized. 

We estimated the over-subtraction by measuring the fluxes 
again in the residual image. Then we re-iterated the global 
subtraction by adjusting the subtracted fluxes and again mea- 
sured the over-subtraction in the next residual image. In each 
iteration, the adjustment of fluxes were small fractions (equal 
to the REALCLEAN gains) of the estimated over-subtraction. 
This slowly brought the residual fluxes of galaxies in the sub- 
tracted image from negative (after the first global subtraction 
following the REALCLEAN cycles) to nearly zero after tens 
of iterations. In the iterations, we closely watched the resid- 
ual fluxes of galaxies to see if they converged to zero. When 
the residual fluxes no longer converged (typically after 40- 
60 iterations), we stopped the iterations. Fluxes subtracted in 
the REALCLEAN stage and in this minimization stage were 
added to make the total IRAC fluxes of the WIRCam K s de- 
tected galaxies. 

For each IRAC band, there are two epochs of observations 
covering different areas in the GOODS-N but with some over- 
lap. All the procedures described above, including the PSF 
construction, were carried out independently on each single- 
epoch image. This is because the two epochs have different 
orientations of pixels and PSFs. For sources in the overlap 
region, we adopted the error weighted fluxes from the two 
epochs. 

6.1.5. Flux and Error Measurements 

As mentioned, we used 4" -6" diameter apertures to mea- 
sure fluxes, or in other words, to evaluate whether the model 
best removes fluxes of galaxies. Such small aperture sizes 
nicely cover the cores of the IRAC PSFs and minimize in- 
fluences from nearby sources. We also derived aperture cor- 
rections from the convolution kernels. However, the proper 
correction is a function of the sizes of the galaxies, which are 
not exactly point sources under the IRAC resolution. This bi- 
ases the measured fluxes, especially when the apertures are 
small. Our small apertures also miss the information in the 
outer wings of the PSFs. To account for these, we repeated 
the REALCLEAN procedures with 50% larger apertures and 
compared the results with those from small-aperture REAL- 
CLEANs. Such aperture sizes are comparable to the sizes 
of the first diffraction rings in the IRAC PSF and are suffi- 
ciently large for the majority of galaxies in the GOODS-N. 
We applied a constant correction factor (2%-6%) to all the 
small-aperture fluxes, derived from the large-aperture REAL- 
CLEAN on high S/N and isolated sources. 

We estimated flux errors in a similar way to those for the 
WIRCam fluxes described in § We convolved the residual 
IRAC images with circular uniform kernels of amplitude 1 .0 
and measured the rms fluctuation locally. This way, the mea- 
sured noise included not only photon noise, confusion noise, 
and correlated noise, but also rough estimates of the uncer- 
tainties caused by imperfect object removal in the entire RE- 
ALCLEAN and flux minimization processes. IRAC fluxes 
and errors measured in this way are included in the Ks se- 
lected source catalog in this data release (TableO. 

6.2. Quality and Performance 
6.2. 1 . Residual Images 



To evaluate the performance of the above REALCLEAN- 
based procedures, we first visually inspected the residual im- 
ages after subtracting the best models. Figure [8] shows the 
HDF-N region of the IRAC images before and after the sub- 
traction. Overall, the subtraction of IRAC objects with our 
REALCLEAN method appears reasonably good in all IRAC 
channels. However, minor residual effects can be seen. There 
are faint halos around many objects in all channels after the 
subtraction, likely a consequence of the fact that our small 
flux apertures did not fully account for fluxes from slightly 
extended objects. On the other hand, on compact objects 
there are residual patterns around the cores of the objects. 
In the first two IRAC channels, this indicates a ~ 5% uncer- 
tainty in the adopted IRAC PSF. Because this occurs across 
the entire field, this is unlikely to be an effect of PSF varia- 
tion or position-dependent mis-registration. This is the most 
severe in the first two IRAC channels, suggesting that under- 
sampling of the PSF by the large IRAC pixels may have 
played a role. However, the exact reason remains unclear to 
us. The 5.8 /im channel has the cleanest subtraction of IRAC 
objects. It gets slightly worse again in the 8.0 /im channel, 
likely because the 8.0 /im morphology becomes sufficiently 
different than that in the /^-band. 

6.2.2. Monte Carlo Simulations 

The above visual inspection of the residual images quali- 
tatively shows that overall our REALCLEAN approach can 
account for the IRAC fluxes of X^-band detected objects and 
that there are no catastrophic failures in subtracting off the 
IRAC objects. However, this does not guarantee that fluxes 
are correctly assigned to objects when multiple objects are 
blended by the IRAC PSF. To test this, we carried out simple 
Monte Carlo simulations. 

We simulated objects of high S/N (~ 20) blended by nearby 
objects. The flux ratios between the blending bright sources 
and the targets are from 1.0 to 100. The separations between 
the targets and the nearby sources are from l."3 to 4", with 
various position angles. All objects are point-like and have 
flat spectra. PSFs derived from the two epochs of the real 
IRAC images and the real WIRCam image were used to cre- 
ate the associated simulated images. All the images were 
processed identically to the above REALCLEAN and catalog 
making procedures. 

We shown the results of the simulations in Figure [9] A few 
trends can be observed and are indeed as expected. REAL- 
CLEAN performs the best when the fluxes ratios are < 10 and 
separations are > 2". At very small separations of l."3, RE- 
ALCLEAN can also correctly recover the fluxes if the nearby 
objects are less than 2x brighter then the targets. Among the 
~ 15 000 r eal sources detected at 3.6 /im in our final cata- 
log (§ 16. 3b . 90% of their brightest neighbors fall in the above 
ranges of flux ratio and angular separation. 

There no data points in regions of large flux ratios and small 
separations (including flux ratios of greater than 31.6, which 
are not shown). These are the cases in which even SExtrac- 
tor cannot resolve the blended sources in the WIRCam im- 
ages. REALCLEAN cannot determine the IRAC fluxes if the 
sources are not resolved by SExtractor in the WIRCam im- 
ages. At the boundary when the SExtractor can barely resolve 
the blended pairs, REALCLEAN either over-estimates (large 
separation) or under-estimates (small separation) the fluxes of 
the targets. The amount of over and under-estimate can be up 
to 40%. Here we are pushing the limits of both REALCLEAN 
and SExtractor. In our catalog, the fractions of sources with 
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FIG. 8. — Results of the REALCLEAN subtraction and residual flux minimization. Images to the left are the original IRAC images in the HDF-N. Images to 
the right are the residual images after subtracting all sources with minimal residual fluxes. From top are the IRAC images at 3.6, 4.5, 5.8 and 8.0 fim. Each 
residual image is shown with a brightness scale identical to its original image. All images are shown with inverted scales. 



< 4" neighbors that are > 17.8 x and > 31.6x brighter are 
5.7% and 3.6%, respectively. 

The above results imply good REALCLEAN fluxes for the 
majority of sources in our catalog. A small number of sources 
may have IRAC fluxes affected by very bright neighbors. The 
effect is generally not ^> 40% before SExtractor becomes un- 
able to detect the faint sources. The caveat here is that the 
above simulations are very simple. One may conduct much 
larger simulations with variable object morphology, S/N level, 
and number of blended sources. We do not think these are par- 
ticularly usefu l given the large uncertainty in the IRAC PSF 
discussed in § 16.2.11 Realistic simulations should take into 
account such uncertainties, folding in how the PSF is under- 
sampled by the IRAC pixels and how the observations were 
carried out and reduced to better sample the images (e.g., 
dither and drizzle), which is out of the scope of this paper. 
Without good knowledge about the PSF, going further beyond 
the above simulations does not likely help on understanding 
the performance of REALCLEAN. In what follows, we will 
rely on simple consistency checks within the available data, 



paying special attention to blended galaxies. 

6.2.3. Overlap Region Between Epochs 1 and 2 

The GOODS-N IRAC observations consist of data from 
two epochs. Both the PSFs and the pixels in these two epochs 
have different orientations. The two epochs of imaging over- 
lap in approximately 35 arcmin 2 of area at the center of the 
GOODS-N. Since we ran our REALCLEAN procedures on 
the two epochs separately, we compared our results in the 
overlap region to see if they are consistent. Figure [TOl shows 
the ratios of the fluxes derived from epochs 1 and 2. The 
dashed curves in each plot show the dispersion of the flux ra- 
tios, and the solid curves s how th e expected dispersion based 
on flux errors estimated in § 16.1.51 The fact that the two curves 
are very close to each other is a strong indication that our er- 
ror estimate is realistic and that the uncertainty of the fluxes 
derived from our REALCLEAN procedure is limited by the 
noise in the data. 

We examined the IRAC sources that have nearby A^-band 
detected objects with separations of l"-2" (diamonds and 2"- 
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FIG. 9. — Simple Monte Carlo simulations of REALCLEAN on two blended objects. Four large panels show results from the four IRAC channels. The 
simulations include various flux ratios and angular separations for the blended pairs. Results of different flux ratios (FR) are shown in the sub-panels. Vertical 
axes show the output-to-input flux ratios of the fainter objects between the blended pairs. Error bars show the dispersion among various realizations. Note that 
the 8.0 fim panel has a different vertical scale. The lack of data points in some regions of the plots means that even the SExtractor cannot separate the blended 
sources in the WIRCam images. REALCLEAN cannot deal with such cases since it is based on the SExtractor WIRCam catalogs. 



3" (triangles). In the l"-2" category, the distributions of the 
flux ratios are not significantly different from the main sam- 
ple, except that at 8.0 fim. At 8.0 /im the 2"— 3" category 
has a distribution similar to the main sample. These results 
suggest that our procedures produce consistent results against 
slight changes in the PSF for blended sources with separa- 
tions greater than 1" at 3.6, 4.5, and 5.8 /mi and 2" at 8.0 /an. 
In the Ks catalog for sources with S/N > 5, the median sep- 
aration between galaxies and their nearest neighbors is 4."1, 
and 99.5% of sources do not have a nearest neighbor within 
1". This suggests that nearly all of the fluxes we derived for 
the first three IRAC channels are reliable. On the other hand, 
among the ~ 5000 8.0 /mi detected sources, approximately 
10% of the sources have nearest neighbors at < 2". Their 
fluxes are more problematic, if Figure[l0]is indicative. How- 



ever, the percentage would be lower if we only considered 
sources whose nearest neighbors were significantly brighter 
than themselves. One can imagine that if sources with similar 
fluxes are blended, REALCLEAN naturally assigns similar 
fluxes to both if the gain is sufficiently small (which is the 
case for our 8.0 /tm REALCLEAN). On the other hand, if 
sources with different fluxes are blended, the procedure might 
fail, especially on the fainter member of the pair. These are 
also indicated in our Monte Carlo simulations. 

6.2.4. Comparison with SED Models 

The key goal of our analyses is to obtain accurate colors 
for the WIRCam and IRAC sources. We selected sources 
that have redshifts in the spectroscopic survey of B08 and 
compared their colors with model SEDs. Here we used two 
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FIG. 10. — Comparison of the IRAC fluxes derived from the two epochs. Dots show the full sample. Diamonds and triangles show galaxies whose nearest 
Kg neighbors have angular separations of 1 "-2" and 2"-3 " , res pectively. Dashed curves show the dispersion of the full-sample data. Solid curves show the 
expected dispersion based on flux uncertainties derived in § 16.1.51 The y-axes have logarithmic scales. 



SEDs with dust and polycyclic a romatic hydrocarbon (PAH) 
features from ISilva et alj dl998l) . M 82 (a dusty starburst) 
and M 100 (a spiral ), and a dustless starburst template from 
iKinnev et al.l (Q~996) that has a blue stellar continuum and 
strong recombination li nes. We also reddene d the M 82 SED 
with Ay = 2.0 using the ICalzetti et al.l (120001) extinction law. 
These four SEDs represent a broad range of galaxy colors. 
Figure QT| shows the comparison between the data and the 
model SEDs. As in Figure [10] galaxies with nearest neigh- 
bors of l"-2" separations are shown with diamonds. 

Overall the derived IRAC and Ks colors match the SED 
models excellently. At z ~ 2.2 and 3.4 in the Fj^/Fks ratio 
and z ~ 3 in the F4.5 /F3.6 ratio, there are substantial numbers 
of observed colors significantly bluer than the starburst model. 
This is true for both isolated galaxies and galaxies with close 
neighbors. We believe these blue "outliers" are not artifacts 
of our REALCLEAN procedures but rather objects with emis- 
sion lines even stronger than those in the starburst template. 
The lines are likely Ha and [S II] (entering Ks at z ~ 2.2), H/3 
and [O III] (entering K s at z ~ 3 .4), and [S III] (entering the 3 .6 



fjm band at z ~ 3), all present in the starburst template. We 
attempted to verify this with our spectroscopic data. Unfortu- 
nately, we only have optical spectra on th ree of these galaxies , 
as most of the redshifts are adopted from lReddv et al.l (2006). 
The three sources have standard Lyman break galaxy spectra 
with strong C IV absorption. In one case there is also Lya 
emission. In the other cases we do not cover Lya. Galaxies 
with strong rest-frame optical (corresponding to Ks and IRAC 
bands at these redshifts) emission lines have absorption lines 
in the rest-frame ultraviolet (corresponding to optical at these 
redshifts). Given their rest-frame ultraviolet continua, these 
three galaxies are inevitably strong emission-line galaxies in 
the rest-frame optical. This is consistent with strong emission 
lines causing the blue observed colors at Ks and IRAC. 

There are a few unexplained blue outliers at z ~ 1 .4—2 in 
Fi.s/Fks and at z ~ 2 in F4.5/F3.6. They might be fluxes mis- 
takenly assigned by our REALCLEAN procedure. However, 
if this is the case, it would be very surprising if such outliers 
are not observed at all at z < 1.5, where most of the spec- 
troscopic samples lie. An alternative explanation is that they 
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are galaxies with misidentified redshifts. We have spectra for 
these galaxies, except fo r two whose redshifts were adopted 
from lReddv et al.l (1200 6). We inspected the available spectra 
and did not see obvious signs of misidentification. Thus, the 
nature of these blue outliers remains unclear. Nevertheless, 
we emphasize that the observed IRAC colors of the majority 
of the galaxies match excellently with the models at nearly 
all redshifts. Combining this with our observations in § 16.2.31 
we are confident that our REALCLEAN procedures correctly 
assigned fluxes to nearly all of the IRAC sources. 

6.2.5. Calibration 

The last performance test is a comparis on with Spitzer ab- 
solute flux calibrators in Reac h et alj ( 120051) . Figure [T2lshows 
the measured colors of the spectroscopically identified stars 
brighter than Ks = 19 and the colors of the IRAC primary ab- 
solute calibrators in Reach et al. Among the calibrators, A 
dwarfs are shown with filled squares and K giants are shown 
with open diamonds. Reach et al. found a displacement be- 
tween the calibration derived from A and K stars, especially in 
the first two IRAC channels. The reason for this is unclear and 
the authors decided to adopt an absolute calibration based on 
only A stars. In Figure [12] we see a gradual drift of the IRAC 
colors for our data with respect to the IRAC calibrators, from 



blue to red in the F^/Fks color. However, since the IRAC 
absolute calibration is based on A stars, it is more appropriate 
for us to compare A stars only. At the very blue end where we 
find three of the four A-type IRAC calibrators, our observed 
colors for the stars agree excellently with the IRAC calibra- 
tors. The largest offset is observed in the F4.5/F3 6 color (top 
panel), where the stars in the left appear to be ~ 0.05 magni- 
tude bluer than the mean color of the A-type calibrators (mea- 
sured from stars with.x-axis values < -0.9). We thus conclude 
that our colors are well calibrated, with systematic errors no 
larger than ~ 5%. 

6.2.6. Known Issues 

Despite all the tests above, the IRAC fluxes derived by our 
REALCLEAN procedures are not perfect. The PSFs we used 
for REALCLEAN can only effectively remove fluxes within 
~ 4x the FWHM. There are about 10 or so very bright stars 
with PSF wings, which are much broader than this. Objects 
close to these bright stars, especially those on the diffrac- 
tion spikes and crosstalk features, have systematically higher 
IRAC fluxes. A more subtle effect is seen in a few very 
crowded regions where more than three objects heavily over- 
lap. In such configurations objects that are close to the centers 
of the groups of objects can have systematically higher fluxes. 
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FIG. 12. — IRAC colors of spectroscopically identified stars (dots) and 
IRAC primary absolute calibrators. A-type calibrators are shown as filled 
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Visual inspections identify only a handful of such examples in 
each IRAC band. 

6.3. Combined NIR and IRAC Catalog 

We combined the IRAC fluxes derived with the REAL- 
CLEAN procedures with the K$ source catalog to form the 
final multiband catalog. The catalog is presented in Table [3] 
The total number of sources is identical to the SExtractor Ks 
catalog, but not every Ks source contained a detection in the 
IRAC bands. As is the case fo r the K s fluxes, our IRAC flux 
errors are quite realistic (see § 16.2.31 and Figure fTTb . and we 
consider the IRAC sources with S/N < 1 to be undetected. Ta- 
ble 2] summarizes the number of detections and upper limits 
in the multiband catalog. 

We have searched for missing sources in the REAL- 
CLEANed IRAC images. There may exist sources that are 
so red that they are detected in the IRAC channels but not in 
the Ks image. These sources are mostly too faint to be de- 
tected in the WIRCam Ks image. A small fraction of them are 
blended with nearby sources in the Ks image and not extracted 
by SExtractor. To search for such sources, we combined the 
REALCLEANed (those in the right hand side of Figure© 3.6 
and 4.5 /im IRAC images and masked the Ks detected objects. 
This naturally excludes sources that are too close to bright 
objects. We visually inspected the REALCLEANed images 
and look for rermaining sources there. We excluded sources 
that are likely cosmic rays in the IRAC images, or other kinds 
of artifacts. We also examined the 5.6 and 8.0 /im REAL- 
CLEANed images, and we found that all the sources seen in 
these two images can be extracted from the 3.6+4.5 /im RE- 
ALCLEANed images. Thus, the positions of such sources 
were measured in the 3.6+4.5 /im image and their fluxes were 
measured in all the individual IRAC images. The A^-band 
fluxes and flux errors are measured in the same way as de- 



scribed in § The positions, Ks fluxes, and IRAC fluxes of 
a total of 358 such sources are listed in Table [5] All of them 
have detections in at least two IRAC bands, and 1 12 of them 
are detected in all four bands. As can be seen in Figure [8] 
source extractions in the REALCLEANed images are difficult 
and can be severel y affec ted by the imperfect source subtrac- 
tion described in § 16.2.11 We warn the readers that Table |5]is 
very likely an incomplete list of X^-faint IRAC sources. 

6.4. Alternative Color Measurements with 
Cross- Convolution 

In the previous subsections, we described our REAL- 
CLEAN method of obtaining robust color measurements, for 
overcoming the issues of the very different PSFs and the 
blending of sources. Here we describe an alternative method 
that can better handle the former issue but is more vulner- 
able to the latter. We started with the background sub- 
tracted and registered images used for the REALCLEAN pro- 
cedure. For each single-epoch IRAC image, we convolved 
the WIRCam K s image with the IRAC PSF and the IRAC 
image with the WIRCam PSF. We refer to this procedure as 
cross-convolution (XCONV). We then used the 2-image mode 
of SExtractor to detect sources in the unconvolved WIRCam 
image and to measure the WIRCam and IRAC fluxes in the 
XCONVed images at the locations of the detected WIRCam 
sources. The PSFs in the XCONVed images have FWHMs of 
slightly larger than 2". After experimenting with various pho- 
tometry apertures, we found that a 3" aperture gives the best 
compromise for S/N, inclusion of most fluxes from a wide 
range of source structure, and blending of sources. Here we 
focus on the /^-to-IRAC colors measured this way, and we do 
not attempt to measure total fluxes. 

The advantage of the XCONV method is clear: the PSFs 
are ideally matched and deconvolution of the IRAC PSFs 
is not required (unlike in the REALCLEAN method). The 
disadvantage is the severely degraded Ks image resolution, 
which prevents the extraction of information on both blended 
sources and faint sources. Because of these, we consider 
the XCONV colors of relatively bright and isolated sources 
to be more robust than the REALCLEAN colors. However, 
as shown in Figure [13] the XCONV colors and the REAL- 
CLEAN colors turn out to be remarkably consistent. The 
systematic offsets in Figure Qj] between these two colors for 
- 1800 F Ks > 10 /iJy objects are 0.002, -0.024, 0.031, and 
0.030 (from 3.6 to 8.0 /im). We examined the small num- 
ber of objects in the figure that have offsets of > 0.05 at 
Ffc s > 10 /iJy. They are either close to bright objects (a weak- 
ness of the XCONV method) or are in regions with multi- 
ple intermediate brightness objects (a weak ness of the RE- 
ALCLEAN method, as mentioned in § 16.2.61 ). 

The excellent consistency between the XCONV colors 
and the REALCLEAN colors further shows that our RE- 
ALCLEAN method is robust. We thus consider the REAL- 
CLEAN catalog as the primary product of this work, and we 
recommend it for most applications. On the other hand, since 
the XCONV colors have unique advantages in their own right, 
we also provide the XCONV colors in Table [6] We remind 
the reader that the XCONV colors are less reliable on faint or 
blended sources. Such colors should be adopted with caution. 

7. PROPERTIES OF THE Ks AND IRAC CATALOG 

Here we provide general descriptions of the properties of 
our Ks and IRAC REALCLEAN catalogs. In subsequent pa- 
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published in its entirety in the electronic edition of the Astrophysical Journal. A portion 


is shown here for g 
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regarding its form and content. If the measured fluxes are less than the flux errors (S/N < 1), zeros are assigned to the fluxes but the errors are 
given for upper limits. Zero flux errors mean that the sources do not have measurements (i.e., the Ks sources are not in the IRAC regions). 
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FIG. 13. — Comparison between the colors derived with the XCONV and REALCLEAN methods. All sources in both catalogs are plotted. 



pers we will present scientific analyses on high-redshift galax- 
ies based on this catalog and other multiwavelength data. 

7.1. Overlap with Other Catalogs 

The GOODS-N has extremely rich multiwavelength data 
and several catalogs have been published of sources in this re- 



gion. Here we summarize the overlap between these catalogs 
and our Ks selected catalog. 

In the GOODS-N and its flanking fields, the two most im- 
portant opti cal observations are the HST ACS imaging of the 
GOODS-N (iGiavalisco et al.l |20Q4j) and the Subaru Suprime- 
Cam imaging of the ~ 0.45 deg field around the GOODS-N 
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TABLE 4 

Number of Sources and Detection Limits in the Multi-Band 
Catalog 



Waveband 


Total Number 


Detections 


Upper Limits 


Depth (//Jy) 


K S 


94951 


75187 


1784 


0.18" 


3.6 fjm 


16950 


15018 


731 


0.11 


4.5 /im 


16631 


14015 


1000 


0.12 


5.8 /im 


16870 


8481 


4551 


0.42 


8.0 ^im 


15184 


7245 


4118 


0.42 



FIG. 14. — ^5 and IRAC color-magnitude diagrams of IRAC selected sources. Dots show all sources selected at each corresponding IRAC band at > 3cr. Gray 
squares show spectroscopically identified stars. In each panel, solid lines show the median 3<r detection limits of the two wavebands in consideration, adopted 
from Table|4] 

IRAC counterparts (489 sources). The later case may be an 
indication of t he existence of spurious sources in our K$ cata- 
log (see § 15.2b. 

We can also reverse the direction of the counterpart search. 
Among the 39432 ACS sources, 16982 (43.1%) have at least 
one Ks counterpart within 0."8. The number decreases to 
14853 (37.7%) for 0."4 search radii. These show not only 
that the Ks depth is insufficient for detecting most of the ACS 
sources, but also that the Ks resolution is too low for separat- 
in g some close pairs of galaxies in the ACS image. 

Capa k et aH (120041) released a Subaru SuprimeCam catalog 
of the GOODS-N region covering an area of 713 arcmin 2 with 
48858 Rc+Z 1 selected sources. In the same region, there are 
51649 Ks selected sources with S/N > 3, among which 73.2% 
have c ounterparts in the S ubaru catalog within l."2. The cat- 
alog of ICapak et all (2004) adopted a tighter S/N cut than our 
Ks catalog. This partially explains why the fraction with op- 
tical counterparts here is much lower than in the ACS case, 
despite the fact that the Subaru images are very deep. The 
large scattering halos around bright stars in the Subaru im- 
age create "holes" in the catalog, and Ks sources do not have 
optical counterparts in such holes. 

The X-ray Chandra Deep Field-North (CDF-N) is fully 
cove red by our ^ imag ing. The Chandra 2 Ms imag- 
ing (lAlexander et al.1 120031) covers an area of approximately 
420 arcmin 2 with an uneven sensitivity distribution. There 
are 503 X-ray sources in the 2 Ms catalog, 493 (98.0%) of 
which have Ks source counterparts within 2." 5 search radii. 
Among the 10 A^-faint sources, four are close to the edge of 



NOTE. — The total number column gives the number of IRAC flux mea- 
surements at the locations of the Kg objects. They are mainly determined by 
the area coverage of the IRAC images. Detections are measurements with 
> 3cr fluxes. Upper limits are measurements with < \a fluxes. Depths are 
median \u errors among sources detected at 3<t for the entire maps. 
a If we limit the Kg sample to the IRAC area, then the median ltr error 
decreases to 0.12 /ijy. 



(ICapak etalJ 120041) . In the 184 arcmin 2 HST ACS region, 
there are 16326 K s selected sources with S/N > 3 in our 
catalog and 39432 zf (F850LP) selected ACS sources in the 
GOODS v2.0 catalog. With 0."4 search radii (0.5 x K s seeing 
FWHM), 14291 (87.5%) K s sources have at least one ACS 
counterpart. The number increases to 14883 (91.2%) when 
the search radii are 0"8. It increases much more slowly be- 
yond this. Among the 1443 optically faint K s sources, 547 of 
them are detected at and 3.6 /im at > 3<j. These are extremely 
unlikely spurious sources and appear to be truly optically faint 
NIR objects. The rest of the optically faint Ks sources are ei- 
ther outside the IRAC region (407 sources) or do not have 
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TABLE 5 

IRAC Sources without K s Counterparts 



RA. 


Dec. 


Fk, 


OKs 


F3.6 




F4.5 0-4.5 


Fs.s 




Fg.o 


fg.O 


(J2000.0) 


(J2000.0) 










(MJy) 











189.049062 
189.086443 
189.049452 
189.209527 
189.215522 
189.241636 
189.246272 
189.087076 
188.998232 
189.018589 
189.045539 
188.981101 



62.133903 
62.134109 
62.138386 
62.140676 
62.144264 
62.145789 
62.146429 
62.150220 
62.150675 
62.153052 
62.154243 
62.154402 



0.227 
0.000 
0.000 
0.000 
0.000 
0.250 
0.217 
0.000 
0.000 
0.192 
0.000 
0.254 



0.091 
0.108 
0.074 
0.114 
0.094 
0.082 
0.105 
0.087 
0.095 
0.101 
0.093 
0.081 



0.945 
0.400 
0.377 
0.908 
0.638 
0.423 
0.468 
0.756 
0.269 
0.345 
0.572 
1.280 



0.128 
0.127 
0.078 
0.181 
0.212 
0.169 
0.151 
0.232 
0.117 
0.112 
0.129 
0.198 



0.894 
0.480 
0.110 
1.037 
1.049 
0.589 
0.623 
0.366 
0.581 
0.277 
0.511 
2.046 



0.160 
0.151 
0.080 
0.191 
0.190 
0.118 
0.122 
0.211 
0.159 
0.129 
0.146 
0.210 



1.597 
0.000 
1.870 
1.419 
0.000 
0.590 
0.377 
0.696 
0.000 
0.000 
1.249 
3.184 



0.625 
0.247 
0.560 
0.731 
0.284 
0.402 
0.324 
0.300 
0.337 
0.496 
0.549 
0.922 



1.049 
0.000 
0.741 
1.270 
1.507 
0.000 
0.000 
0.000 
0.426 
0.769 
0.000 
2.433 



0.385 
0.614 
0.377 
0.541 
0.590 
0.282 
0.258 
0.425 
0.335 
0.557 
0.438 
1.231 



NOTE. — Table[3]is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is 
shown here for guidance regarding its form and content. If the measured fluxes are less than the flux errors (S/N 
< 1), then zeros are assigned to the fluxes, but the errors are given for upper limits. 



TABLE 6 

K5-T0-IRAC Colors Derived with Cross-Convolution 



ID 


Kg -3.6/im 


K s -4.5p.m 


Ks-5.8/im 


Ks-8.0pm 


54739 


-0.839 


-1.434 


-1.807 


-2.430 


54740 


-0.135 


-0.628 


-0.576 


-0.884 


54741 


99.000 


99.000 


99.000 


99.000 


54742 


-0.579 


-0.928 


-1.283 


-1.879 


54743 


-0.517 


-0.862 


-1.095 


-1.537 


54744 


99.000 


99.000 


99.000 


99.000 


54745 


0.154 


-0.067 


0.166 


0.174 


54746 


99.000 


99.000 


99.000 


99.000 


54747 


99.000 


99.000 


99.000 


99.000 


54748 


0.365 


0.190 


-0.273 


-0.473 


54749 


-0.659 


-0.817 


99.000 


99.000 


54750 


0.283 


-0.006 


-0.249 


-0.316 



NOTE. — Table [6] is published in its entirety in the electronic 
edition of the Astrophysical Journal. A portion (identical to that for 
Table[3j is shown here for guidance regarding its form and content. 
Source ID and the total number of sources are both identical to that 
in Table [3] All colors are in the AB magnitude system. Sources 
with no measurements are assigned 99.0. 



relatively bright objects. The brighter Ks objects may have af- 
fected the detection of the X-ray counterparts at Ks, or them- 
selves may be the correct counterparts. To check whether the 
rest 6 are spurious X-ray sources, we look for their counter- 
parts at other wavelengths. Two A^-faint X-ray source ( source 
#151 and #183 in the catalog of lAlexander et al.ll2003h have 
faint IRAC counterparts and are extremely faint in the optical. 
They are even close to two submilli meter sources detec ted at 
1.1 mm (AzGNIO amd AzGN12 in lPerera et all 12008). An- 
other (#88) has a faint optical counterpart but is outside the 
IRAC region. These three are therefore genuinely A^-faint X- 
ray sources. Two A^-faint X-ray sources (#248 and #371) are 
undetected in the Subaru optical, Spitzer 3.6 to 24 /an, and 
VLA 1.4 GHz images. Even in the X-ray images they ap- 
pear to be quite faint. They are more likely spurious X-rayt 
sources. Finally, one X-ray source (#427) is not detected in 
any of the optical to radio images we have, but its detection 
appears quite convincing in the hard X-ray band. Its nature is 
undetermined. 

The Spitzer GOODS Legacy Program also imaged the 
GOODS-N fiel d with Multi-Ban d Imaging Photometer for 
Spitzer (MIPS. iRieke et all 120041) at 24 /im. The GOODS 
Data Release 1+ provided a 24 /mi catalog that is highly com- 
plete at 80 /iJy, containing 1 199 sources. With search radii of 



2."5, 1 182 MIPS sources have counterparts in our Ks catalog. 
The rest are all in low S/N regions in the MIPS images and 
appear to be spurious. Thus, our Ks image seems to detect all 
> 80 /iJy 24 /im sources, bearing in mind the possibility that 
some sources could be misidentified by the large search radii. 

Finally, at the long wavelength side, there exist dee p 
VLA 1.4 GHz imagin g in the GOODS-N (lRichardsll2000l) . 
lBiggs&Ivisonl (2006) re-reduced the VLA data and provided 
improved images and catalog. The catalog of Biggs & Ivison 
contains 537 > 30 /iJy (~ 5er) 1.4 GHz sources within a ra- 
dius of ~ 24' centered at the HDF-N with a non-uniform sen- 
sitivity distribution. A total of 450 radio sources are within 
the 30' region of our K$ image, 439 (97.5%) of which have 
Ks counterparts within 1 "5 of search radii. The 1 1 ^-faint 
radio sources do not appear to be spurious in the radio im- 
age, although some of them are very close to the 5cr limit at 
1.4 GHz. They are all in the outer shallower region of the Ks 
image. We note that, in the outer part of the image, the Ks sen- 
sitivity decreases dramatically with exposure time (Figure[2]), 
while the radio sensitivity decreases slowly with the primary 
beam of the VLA. We therefore conclude that the central part 
of our Ks image is sufficiently deep for detecting nearly all of 
the 30 /iJy radio sources. 

7.2. Color-Magnitude Diagrams 

In Figure [14] we present color-magnitude diagrams for 
sources selected in each of the four IRAC bands. The most 
apparent feature is the color bimodality in the F%g/Fks and 
F4.5/F3.6 colors. There are two distinct, nearly horizontal 
branches. The blue branches contain most of the stars and 
low-redshift galaxies. The redshift divisions in F-ss/Fks and 
F4.5/F3.6 are z ~ 0.7 and z ~ 1.2, respectively. At these two 
redshifts, galaxies quickly move across the two branches, be- 
cause of the rest-frame 1.6 /im SED bump in galaxies. This 
can also be seen in Figure [TT] which is based on our spec- 
troscopic redshifts. The division at z ~ 1.2 in the F4.5/F3.6 
color is much more apparent in Figure [14] because Figure [Til 
suffers a bit from the spectroscopic redshift desert at z > 1 .4. 
The bimodality nearly disappears in the F5.8/F4.5 color and 
entirely disappears in the Fg.o/Fs.g color. The main reason is 
that strong PAH emission from low-redshift starbursts signif- 
icantly reddens the colors. A minor reason is the relatively 
lower S/N at these longer IRAC bands. 

In Figure [15] we show the color-magnitude diagrams for 
the X-ray, 24 /im, and 1.4 GHz selected sources discussed 
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Fig. 15.- 
10 42 erg s" 



- Ks and IRAC color-magnitude diagrams of X-ray, 24 fim, and 1 
are enclosed by squares. Solid lines in each panel show median 3 



4 GHz selected sources. In panel (A) and (B), X-ray sources more luminous than 
a detection limits of the two wavebands in consideration. 



in the previous subsections. The F^^/Fks colors of all three 
populations lie in a well-defined area of the color-magnitude 
space. The distributions are determined by the sensitivities at 
X-ray, 24 /im and 1 .4 GHz, plus loose correlations between 
luminosities at these wavebands (tracing star formation and 
black hole accretion) and the rest-frame 1—4 /im luminosities 
(tracing stellar mass). At bright magnitudes the sources are at 
low redshifts (the blue branch), while at the fainter fluxes the 
sources are at z > 0.7 and the sources move to the red branch. 
The trend is the weakest for X-ray AGNs (open squares in 
Figure \l5k), but it is still observable. 

The distributions in the F4.5/F3.6 and F5.8/F4.5 colors of 
these three populations are more or less similar, and we do not 
show them here. However, the Fs.o /F5.8 color-magnitude dia- 
grams become significantly different. A group of red Fg .0/F5.8 
sources appear at the bright end. These colors are produced 
by strong PAH emission from low-redshift starbursting galax- 
ies. They are also seen in the general 8 /im selected dia- 
gram (Figure [l"4tl). It is interesting to note the difference 
between the low-luminosity X-ray sources (X-ray luminos- 
ity < 10 42 erg s" 1 , dots in Figure [T3b) and the X-ray AGNs 
(> 10 42 erg s" 1 , squares in Figure [T5b). In the Chandra 2 Ms 
sample, none of the X-ray AGNs have red F8.0/F5.8 colors of 
> 0.8, because the strong radiation field of AGNs can destroy 
PAH molecules. 

7.3. IRAC Color-Color Diagrams 

In B08 we presented F4.5/F3.6 vs. F8.0/F5.8 and F5.8/F5.6 vs. 
F8.0/F4.5 color-color diagrams along with our spectroscopic 
sample. With these two IRAC color-color diagrams, we ana- 
lyzed efficiencies and contaminations of selections of z > 1.3 
and z > 1 .6 galaxies as well as AGNs. Since the quality of the 
IRAC photometry is improved here, it is interesting to revisit 
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FIG. 16. — ^4.5/^3. 6 vs. Fg o/Fs.g color-color diagrams of 8.0 fim (> 3<r) 
and X-ray selected sources. In panel (B), X-ray AGNs (> 10 42 erg s~') are 
enclosed with open squares. The dashed lines show the selection of z > 1.3 
galaxies proposed by B08, which is logiF^ s/Fj (,) > 0.096 \og(F% .o/fs. s)- 
The dotted lines show the selection of z > 1.2 galaxies p roposed by 
IDevlin et alj 420091) . The solid lines enclose the region where I Stern et al] 
ism suggested that these colors select most broad-line AGNs. 



these issues. Here we focus on the selections of high-redshift 
galaxies and X-ray AGNs. (B08 also included the most pow- 
erful radio sources, but that sample size is very small.) 
Figure[l6lshows F4.5/F3.6 vs. Fg.o/Fs.g color-color diagrams 
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of 8.0 /im and X-ray selected sources. The color bimodality 
shown in Figure[T4b is clearly seen here. The dashed lines are 
the z > 1.3 selection proposed by B08. We reported in B08 
that this selection can pick up 90% of z > 1.3 galaxies, but 
39% of the galaxies that satisfy this criterion are z < 1.3 con- 
taminants. In the 3a 8.0 /im sample here, there are 2096 spec- 
troscopic redshifts, 352 of which are at z > 1.3. The above 
selection picks up 241 of them, corresponding to a selection 
completeness of 68%, much lower than that in B08. On the 
other hand, the selection only picks up 54 z < 1.3 galaxies, 
corresponding to a selection contamination of 1 8 %, also much 
lower than that in B08. The results do not change much if we 
limit to higher S/N 8 /im sources where the measurements of 
colors are more robust. The inconsistency between the results 
here and our previous results is most likely a systematic dif- 
ference in the calibration of the IRAC fluxes. For example, 
if we make the galaxies here ~ 0.15 magnitude redder in the 
F4.5 /F3.6 color, then the completeness and contamination both 
become similar to those in B08. 

The above comparison shows that using IRAC colors to se- 
lect high-redshift galaxies is extremely sensitive to how galac- 
tic colors are measured. This also emphasizes that simply 
adopting existing IRAC color selection schemes may be prob- 
lematic unless one is very careful about the fine details of the 
color measurements. Bearing this in mind, the selection com- 
pleteness and contamination are then just a matter of fine tun- 
ing the selection criteria at the expense of one another. Such 
fine tuning also requires a relatively large sample of spectro- 
scopic redshifts. In addition to the very likely difference in 
the colors, the basic conclusions about IRAC color selections 
on high-redshift galaxies remain unchanged here and in B08. 

There are also other high-redshift galaxy selections. The 
dotted lines in F igure [TBI show the z > 1 .2 selection adopted 
by iDevlin etail (120091) and iMarsden et al.1 (12009) for studies 
of submillimeter galaxies. We found that all of the above con- 
clusions apply to this selection as well. This selection has ei- 
ther incompleteness at high redshift or contaminations at low 
redshift. 

Figure [T6b shows X-ray sources and X-r ay AGNs. The 
IRAC colors suggested by IStern et all (120051) to select most 
broad-line AGNs are enclosed by solid lines. However, here 
we see that the Stern et al. selection only picks up a frac- 
tion of the X-ray AGNs. By comparing both panels in Fig- 
ure [16] as well as Figures [lTb and d, we see that the same 
region also selects a larger number of high-redshift galaxies. 
T hese galaxie s do no t show up in the relatively shallow data 
of IStern et al.1 (120051 their Figure 1). T his is consistent with 
the finding in B08 and lYun et al] (120081) that high-redshift red- 
dened obj ects occupy simila r color spaces to AGNs. The se- 
lection of IStern et al.1 d2005l) does not pick up the majority of 
AGNs and has substantial galaxy contamination. 

B08 also presented F5 .8/F3.6 vs. F%_q/F^$ color-color dia- 
grams and used them t o analy ze selections of z > 1.6 galax- 
ies and the lLacy et all ((2004) selection of AGNs. As with 
the F4.5/F3.6 vs. F8.0/F5 .8 color-color diagram, B08 concluded 
that the IRAC selection of z > 1.6 galaxies can be very com- 
plete (~ 90%), but it has a large contamination (~ 40%) from 
low-re dshift galaxies. B08 also found that the lLacy et all 
(2004) selection of A GNs is incomplete (though better than 
the lstern et al.l l2~005 selection) and has severe galaxy contam- 
ination. We present the updated F5.8/F3.6 vs. F8.0/F4.5 color- 
color diagrams in Figure \T7\ and we confirm all the conclu- 
sions made by B08. 

To sum up, using IRAC colors to select high-redshift galax- 
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FIG. 17. — F5. 8/^3.6 vs - F$.o/Fa.5 color-color diagrams of 8.0 fim (> 3<r) 
and X-ray selected sources. In panel (B), X-ray AGNs (> 10 42 erg s -1 ) are 
enclosed with open squares. The dashed lines show the selection of z > 
1.6 galaxies proposed by B08, which is log( fsj/ffyi) > 0.09 61og(Fg 0/-F4 .5). 
The solid lines enclose the region where Lacy et al. 12004) suggested that 
these colors select AGNs. 

ies can be either incomplete or suffer from contaminations 
from AGNs or low-redshift galaxies. The selection can be 
fine tuned to balance the completeness and contamination, but 
this requires a good spectroscopic sample and/or extremely 
well measured IRAC colors. We also found that galaxies 
and AGNs occupy a very similar regio n of color s pace. The 
IRAC selections o f AGNs based on the lLacy et"aT] (120041) and 
IStern et ail d2005l) methods always include large numbers of 
galaxy contaminations and have various degrees of incom- 
pleteness. 

7.4. Ks and IRAC Color-Color Diagrams 

The above discussion only includes IRAC colors. In this 
work, we presented ultradeep K$ photometry, which can be 
a powerful tool for studying galactic colors. We show the 
F4.5 /Fk s vs. F4.5 / F3.6 color-color diagrams for 4.5 /im and X- 
ray sources in Figure [18] Here we only adopt the two short- 
est IRAC bands for their higher S/N and for their continuous 
availability in the Spitzer warm phase mission. Again, we 
see that the X-ray AGNs occupy a color space that is nearly 
identical to that of galaxies. Such a diagram is extremely inef- 
fective for separating galaxies and AGNs. On the other hand, 
this color-color diagram is much more structured than the pre- 
vious IRAC-only diagrams, suggesting that these colors can 
provide an excellent redshift diagnostic. 

Figure [19] shows the same diagram as Figure [18] but only 
for the spectroscopic subsample drawn from B08. There is 
a very subtle difference in the distributions of galaxies in 
these two diagrams (4.5 selected vs. spectroscopically 
identified), because the spectroscopic identifications suffer 
from the redshift desert and are somewhat biased to galax- 
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identified sources in the B08 sample. Symbols are color coded with spectro- 
scopic redshifts, including spectroscopically identified stars (gray triangles). 



ies with stro ng emi ssion lines that can have unusual /"i.s/F^ 
colors (see § I6.2.41 i. However, three distinct redshift groups 
can be seen in both diagrams: z < 0.65, z = 0.65-1.2, and 
z > 1.2. In particular, we find that the following selection 
criteria are especially effective in selecting intermediate red- 
shifts of z = 0.65-1.2: 

log(F 4 . 5 /F Ks ) > -0.28 and 
log^s/^.eX-O.l and 
log(F 4 . 5 /F & ) > 1.21og(F 4 . 5 /F 3 . 6 )+0.26. 



This selection is shown as the solid lines in Figure[T9l 

In Figure[l9]a total of 2834 spectroscopic identifications are 
included, among which 1079 are at z = 0.65-1.2. The above 
Ks and IRAC selection picks up 870 galaxies at z = 0.65-1.2, 
corresponding to a selection completeness of 81%. On the 
other hand, the selection also picks up 86 galaxies at z < 0.65, 
26 galaxies at z > 1.2, and 3 stars, corresponding to a con- 
tamination of 11.6%. It is possible to fine tune the selection 
criteria to achieve higher completeness at the expense of in- 
creasing contamination, or vice versa. For example, after we 
moved the selecting region in the diagram towards the lower- 
right by 0. 1 magnitude along both the x and y axes, the se- 
lection completeness increased to 89%, but the contamination 
also increased to 17%. 

Fa.s/Fks and F4.5/F3.6 colors can still be measured in the 
Spitzer warm phase mission and can be used in different ways. 
The above z = 0.65-1 .2 selection picks up galaxies at redshifts 
just lower than th at in the popular BzK selection (z ~ 1 .4-2.5, 
iDaddi et aill2004l) . Galaxies at z > 1.2 and z < 0.65 can also 
be selected from the same diagram. In addition, in a subse- 
quent paper, we will use the same diagram to select z > 2 
massive and dusty galaxies (Wang et al., in preparation). 

Generally speaking, the addition of the Ks-band data in- 
creases the redshift resolution at z < 2. At higher S/N this will 
greatly help IRA C-based photometric redshift studies (e.g., 
iPope et aT]|2006l) . which can be more powerful than simple 
color selections. On the other hand, in all the examples dis- 
cussed in this paper, IRAC colors become much less efficient 
in separating galaxies at z > 2, even after the inclusion of K$. 
For these high-redshift galaxies, optical colors or photomet- 
ric redshifts that com bine optical, NIR, and IRAC bands (e.g., 
Mancini et al. 2009) may be much better. However, on opti- 
cally fai nt objects such as high-redshift subm illimeter galax- 
ies (e.g.. lWang et al.ll2007h lYounger et alJl2009l) . colors from 
deep NIR and IRAC imaging are valuable in studying their 
redshifts. 

8. SUMMARY 

We carried out ultradeep Ks band imaging in the HDFN 
and its flanking field with WIRCam on the CFHT The mosaic 
image has a size of 35' x 35', fully covering the GOODS-N 
field. With nearly 50 hr of integration, we reached a 90% 
completeness limit of 1.38 /iJy (AB = 23.55) in the 0.25 deg 2 
region, and 0.94 /jjy (AB = 23.96) in the 0.06 deg 2 GOODS- 
N IRAC region. The faintest detected objects have Ks fluxes 
of ~ 0.2 //Jy (AB = 25.65). The photometry is uniform across 
the entire field, and the variation is well within 2%. In the 
GOODS-N region the astrometry is accurate to 0"03 for com- 
pact objects with high S/N. The Ks images and the associated 
catalogs are available through the internet. 

We used the Ks image and catalog as priors to measure MIR 
magnitudes in the Spitzer IRAC 3.6 to 8.0 /jm bands with 
our REALCLEAN method. We verified the REALCLEAN 
results using data taken in two different epochs of observa- 
tions, as well as by comparing our colors with the colors of 
stars and with the colors derived from model SEDs of galax- 
ies at various redshifts. The REALCLEAN method appears 
to provide reasonably good photometry even on galaxies that 
have close neighbors within 2". This also provides a rela- 
tively clean way to identify faint IRAC objects that are not 
detected in the K s image. We make the REALCLEAN IRAC 
catalogs for both Kg detected and undetected objects available 
to the community. Additional to the REALCLEAN method, 
we also measured galactic colors by convolving the K s image 
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with the IRAC PSFs and convolving the IRAC images with 
the Ks PSF. The derived colors are consistent with the REAL- 
CLEAN results. 

We found our Ks sample to be extremely powerful in study- 
ing the multiwavelength properties of high-redshift galaxies. 
The Ks image detected most (if not all) of the X-ray, 24 /im, 
and radio sources published in this field. On the other hand, 
only 40% of the HST ACS sources in the GOODS-N are de- 
tected at Ks at the current depth. 

We studied various Ks and IRAC color-magnitude and 
color-color diagrams. Because of the rest-frame 1.6 fim bump 
in the SEDs of galaxies, galaxies show significant color evo- 
lution in the Ks and IRAC bands from z = to z < 2. We 
found that using IRAC colors to select galaxies at certain 
redshifts is extremely sensitive to the color measurements 
and that balancing selection completeness and contamina- 
tion will require careful tuning of the selection criteria with 
a good spectroscopic sample. We discussed an effective se- 
lection of z = 0.65-1 .2 galaxies using the Ks band and the two 



shorter IRAC bands that will remain available in the Spitzer 
warm phase mission. We also confirme d our previous s tudies 
(B08) that the IRA C AGN selections in lLacy et all (12004 and 
IStern et aTl (120051) are either incomplete at the low-luminosity 
end or suffer severely from galaxy contamination. 
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